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Abstract 

We study spin- 1/2 Heisenberg XXX antiferromagnet. The spectrum of the Hamiltonian was found by 
Hans Bethe in 1931, [jjj]. We study the probability of formation of ferromagnetic string in the antiferro- 
magnetic ground state, which we call emptiness formation probability P(n). This is the most fundamental 
correlation function. We prove that for the short strings it can be expressed in terms of the Riemann zeta 
function with odd arguments, logarithm In 2 and rational coefficients. This adds yet another link between 
statistical mechanics and number theory. We have obtained an analytical formula for P(5) for the first time. 
We have also calculated P(n) numerically by the Density Matrix Renormalization Group. The results agree 
quite well with the analytical ones. Furthermore we study asymptotic behavior of P(n) at finite temperature 
by Quantum Monte-Carlo simulation. It also agrees with our previous analytical results. 

1 Introduction 

Recently the considerable progress has been achieved for the exact calculation of correlation functions in 
the spin-1/2 Heisenberg XXZ chain g, |, g, |§ §, @, § |, 0, @. Hans Bethe discovered his Ansatz, 
while diagonalising XXX Hamiltonian. The most important features of exactly solvable models, like two- 
body reducibility of dynamics were first discovered for this model. We believe that the antiferromagnetic 
Heisenberg XXX chain is one of the most fundamental exactly solvable models. Recently we developed a 
new method of evaluation of the multi-integral representation of correlation functions . We study the 

most fundamental correlation function of the model, which we call emptiness formation probability P(n). 
We shall abbreviate it to EFP. It was first introduced in [B| 

n 

P(n) = (GS|JJPi|GS>, (1.1) 

2=1 

where Pj = Sj + \ is the projector on the state with the spin up in the j'-th lattice site. |GS) is the 
antiferromagnetic ground state in the thermodynamic limit constructed by Hulthen P(n) is a probability 
of formation of a ferromagnetic string of length n in |GS). 

The Hamiltonian of the Heisenberg XXX chain is given by 

N 
3=1 



where the coupling constant J is positive for theantiferromagnet. 

In ref. jio) , pl| , we wrote the Hamiltonian (1.2) in terms of the Pauli matrices. It corresponds to J = 4. 
Note that the Hamiltonian (1.2) annihilates the ferromagnetic states f . . . j) (or | J, . . . J.)). 

First cases P(3) and P(4) were calculated by means of the multi-integral representation in |u| [H]]. In 
this paper, we present a new analytic formula for P(5) 



t = !(I) » 1 ,, ((Ilt ™ t(!K(I1 

= 2.011725953 x 10" 6 , (1.3) 
where £(s) is the Riemann zeta function 

oo 

CW = £^> for > 1. (1.4) 

n=l 

P(5) is expressed in terms of In 2, C(3), C(5) and f (7) with rational coefficients. In fact, it was conjectured 
= 

'P(n) is always expressed in terms of logarithm In 2 , Riemann zeta functions C(2fc + 1) with 
odd argument and rational coefficients." 

This means that all values of P(n) are different transcendental numbers Jl4| . 
For comparison, let us list previous results for P(n). 

P(l) = 1 = 0.5, (1.5) 

P(2) = - - -In 2 = 0.102284273, (1.6) 

P(3) = - -In 2+ -C(3) = 0.007624158, (1.7) 
4 8 

P(4) = ^21n2 + ^ C (3)-^ln2. C (3)-|c 2 (3) 

-^C(5) + ^ In 2 • C(5) = 0.000206270. (1.8) 

Let us mention that in contrast with P(l), P(2), P(3) [ which do not contain non-linear terms of the Riemann 
zeta function] the values P(4) and P(5) do contain non-linear terms in Riemann zeta function. It might be 
instructive to express formulae (1.2), (1-5) for P(5) and P(4) in a linear form by introducing the multiple 
zeta values jlf| 

C(ki, ki, . . . , k m ) = 2^ n i kln 2 k2 ■ ■ - n ™ km - (1-9) 

ni>n2>..->n m >0 

The length (or depth) of this multiple zeta value is equal to m and the level (or weight) is equal to fei + 

k-2 + . . . + k m - The result looks as follows: 

1 17^ 11 FiFi «c: 

P($ = i-21n2+-^C(3)-yln2.C(3)-gc(5) + gln2.C(5) 

51 1 5S 

~To c(3 ' 3) + W c(4 ' 2) ' (L10) 

m . >_- t2+ - <( 3 ) _«,, <( 3 ) -^ <(5) 

_-« 3 , 3) + - {(4 , 2) _l^« 4 , 4) _| c(5 , 3) -gc ( a, 2 , 

487037 A/r r , 584425 „ fa „, 596647 _ ,„ „. 42777 A .„ 

H C(5,5) H C(6,4) H C(7,3) H C(8,2). (1.11) 

512 s ^ ' ; 1024 sv ' ; 1024 sv ' ; 1024 sv ' ; y ' 

Here we used the following identities 

C(3) 2 = 8C(3,3) - 3C(4,2), 

C(3)-C(5) = ^C(4,4) + i<(5,3) + ic(6,2), 

C(5) 2 = 22C(5,5) + 10C(6,4) + 10C(7,3), 

C(3)-C(7) = i^C(5,5) + |c(6,4) + ^C(7,3) + ic(8,2). (1.12) 



2 



These identities can be derived using the reccurent relations (4.1), (4.5),..., (4.7) of ref. |lq ]. Let us note 
that the Drinfeld's associator is also related to multiple zeta values in a linear way [jl^ . 

Since, In 2 in formulae above looks somehow isolated, we believe that it seems to be more appropriate to 
express P(n) in terms of the alternating zeta series (the value of polylogarithm at root of unity) 



U") = 



-iy 



n>0 



= -Li s (-1) 



(1.13) 



Here Li s (a;) is the polylogarithm. The alternating zeta series is related to the Riemann zeta function as 
follows 

= i^^Us) (1.14) 

This formula is valid for s / 1. In contrast with the zeta function [which has the pole when s — > 1], the 
alternating zeta has a limit as s — > 1 

Ca(l) = m2 (1.15) 



Using the formulae (1.3),..., (1.8), (1.14) and (1.15), one can get the five first values of P{n) expressed via 
the alternating zeta series 



P(l) 
P(2) 

P(3) 
P(4) 

P(5) 



= g{l-C«(l)}, 

= i{l-4C.(l) + 2C.(3)}, 

= ~{l - 10C.(1) + ^Ca(3) - ^Ca(5) - ^Ca(l) • C.(3) 



<a(3)}, 

1 ( _, nn , 281 . 1355 

_{l-20Ca(l) + — Ca(3) 



+ — Ca(l)-<a(5) ., 



,. Ca(5) + — Ca(7) 
-180Ca(l) ■ C«(3) + ^Ca(l) ■ Ca(5) - ^pCa(l) ■ Ca(7) - ^Ca(3) ■ C«(5) 

+679< a (3) ■ C«(7) - 326< a 2 (3) - ^d(5)}. 



(1.16) 



The values of the Riemann zeta function at odd arguments appear in several places in theoretical physics, 
not to mention in pure mathematics. Transcendental number £(3) first appeared in the expression for 
correlation functions in Takahashi's papers j^j, . He evaluated the second neighbor correlation 



(Si • S l+2 ) = - 



41n2 +i C(3) 



(1.17) 



It was obtained from the 1/U expansion of the ground state energy for the ha lf-fil led Hubbard chain (see 
also a nother derivation in ref. pll]). We remark that the expression for P(3) (1.7) can be extracted from 



(LIZD 



Now let us discuss the asymptotic behavior of P(n) when n is large. At zero temperature we believe 
that P(n) should show a Gaussian decay as n tends to infinity Jl(], [llj]. In order to prove it mathematically 
we have to obtain a general formula for P(n) which has not been achieved yet. On the other hand, we can 
calculate P(n) by numerical means and confirm our analytical results. We have applied the Density Matrix 
Renormalization Group (DMRG) method p3| , ^| and obtained more numerical values of P(n). The result 
is given in the end of section 2. At finite temperature it was shown that P(n) decays exponentially [^o|. 
This time we can employ the Quantum Monte-Carlo (QMC) simulation |^5| to calculate P(n) numerically. 
In section 3 we confirm that P(n) exhibits an exponential decay at finite temperature and confirm again our 
analytical expression. Let us note that this numerical approach was successfully applied to the XX model 
in ref. 01. 



2 P(n) at zero temperature 

The integral representation of P(n) for the XXX chain was obtained in [^) based on the vertex operator 
approach |3| 



dXi f d\ 2 f d\ n Y\{i | I ) n ~ a ( 7r ^ a )" TT sinh7r(Afc - A-,) 
c 2m\i J c 2ni\2 J c 2iri\ n J- A A a sinh7rA a J- A 7r(Afe — Aj — i) 



l<j<k<n 
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The contour C in each integral goes parallel to the real axis with the imaginary part between and — i. 
Below we sketch the evaluation of the integral for P(5). It can be written in the following form 

P(5) = f[ J ^f/ 6 (A lj ...,A 5 )T 6 (A 1 ,...,A 5 ) ) (2.2) 



where 



and 



TT,,, . sinh 7r(A, — At) 

£/ 5 (Al,...,A 5 J = 7T — g — ; , (2.3) 



n S =1 sinh 5 ^A 



r B (Ai,-,A»)= J. 1 '" 1 1 ,\ ' (2.4) 



Taking into account the properties of the functions Ub(Xi, ■ ■ ■ , As), one can reduce the integrand Ts(Ai, . . . , As) 
to the "canonical form", T§ (Ai, . . . , A5) as was done for P(2), P(3) and P(4) in ref. (ll) 

^ (Al <- ' As) ^ S5) + aSt + (a 2 -ao(a 4 -A3) ' ^ 

where the P (0) , P (1) , P (2) are polynomials of the integration variables Ai, . . . , As. The manifest form of these 
polynomials is shown in Appendix. 

Now using the methods developed in |nj we can calculate three integrals that contribute into P(5) 

J o 5) =Il/ fHtMAi,...,As)P (5) (A 1 ,...,As), (2.6) 
j=1 Jo n 



P 2 (j) (Ai, . . . , As) 
A5) (A 2 -A 1 )(A 4 -A 3 )- (2 ' 8) 



The result looks as follows: 



- H (2.9) 
^ " -I-f--^-f« 5 ) + Sc(7,, 
2 12g sw 2 sw lg s,v ; 64 ; 

Summing up these three values we come to our final answer ( |l.3| ). 

Using the same method we can, in principle, get the analytic formula for any P(n). Unfortunately, so 
far we have not succeeded in calculating P(n) for n > 6. On the other hand, one can estimate the numerical 
values of P(n) using DMRG p2| p3|. The method is suitable for studying ground-state properties. We 
followed standard algorithm, which can be found in literature (see ref. p3]). Below we shall outline some 
technical points that are relevant for the simulation precision. We implemented the infinite-system method. 
We have repeated renormalization 200-times. At each renormalization, we kept, at most, 200 relevant 
states for a (new) block; namely, we set m = 200. The density-matrix eigenvalue {w a } of remaining bases 
indicates the statistical weig ht. We found w a > 10~ 10 : That is, we have remained almost all relevant states 
with appreciable statistical weight w a > 10 -10 through numerical renormalizations. In other words, we 
have discarded (disregarded) those states with exceedingly small statistical weight w a < 10~ 10 , which may 
indicate error of the present simulation. 
The obtained data are shown in Table |l| 

Compared with the exact values jl.3| ), jl.5j ),...,(1.8), we see the DMRG data achieve about 3 digits 



accuracy up to P(5). Probably the other data, i.e.,P(6), P(8) also maintain at least 1 or 2 digits accuracy. 
Then combining the exact values up to P(5) and the numerical data in Table 1, we have made a semi-log 
plot in Fig. 1. On the vertical axis we plotted the lnP(n), on horizontal axis we plotted n 2 . From Fig.l we 
clearly see the data fall into a straight line suggesting that the asymptotic form of P(n) is governed by the 
Gaussian form, 

P(n)~a" n2 . (2.12) 
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Table 1: DMRG data for P(n) with uncertainties in the hnal digits. 
~n P(n) by DMRG 



2 


1.0222 x 10" 


— — 

- 1 


• > 

•i 


7.6238 x 10 


-3 


4 


2.060 x 10" 


-4 


5 


2.010 x 10" 


-6 


6 


7.05 x 10" 


-9 


7 


8.85 x 10- 


12 


8 


3.7 x 10" 


15 



We can read off the Gaussian decay rate a from the slope. Our estimate is a — 1.6719 ± 0.0005. It is an 
intriguing open problem to associate this number with a certain analytical expression. Doctor A.G. Abanov 
confirmed Gaussian form of asymptotic expression in the frame of bosonization technique. 



ar-20 




Figure 1: P(n) at zero temperature 



3 P(n) at finite temperature by QMC simulation 

At finite temperature T the emptiness formation probability P(n) is defined by the thermal average 

Tr {e- H/T I\U P >} 



It was shown in ref. 



P(n) Tr {e-v/T} 

that P(n) decays exponentially at finite temperature as n goes to infinity 



(3.1) 



Pin) ~ c (T)e— 



(3.2) 



Here / is the free energy per site of the system and co(T) is a constant prefactor. We shall confirm this 
formula by calculating P(n) numerically by QMC simulation pq|. We adopted the continuous-time algorithm 
|26[ with the cluster-flip update ^?], g^] which is completely free from the Trotter-decomposition error. 
We treated system sizes up to N = 128, and imposed the periodic boundary condition. We performed five- 
million Monte-Carlo steps initiated by 0.5 million steps for reaching thermal equilibrium. P(n) is measured 
over the five-million Monte-Carlo steps. The semi-log plots of the data for several diff eren t temperatures 
are shown in Fig. 2. The straight lines (slopes) are the analytic asymptotic formula (3.2). We assumed 
co(T) = 1 for simplicity. Actually co(T) deviates from unity for low temper atur e. Anyway we observe that as 
n — > oo, EFP P(n) decays exponentially according to the asymptotic form (3.2). In particular, at sufficiently 
high temperatures, even for small n, EFP P(n) is well fitted by (3.2). In contrast, at low tempe ratu re , 
when n becomes small P(n) may reflect the Gaussian decay at zero temperature and deviates from (3.2). It 
agrees with the physical picture. Since, zero temperature is a critical point, we expect qualitative change of 
asymptotic of correlation functions. 

Finally let us make a comment how we can evaluate the free energy per site / from the point of view of 
the Bethe Ansatz. There are now three different integral equations which determine the free energy /; 
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1. Thermodynamic Bethe Ansatz (TBA) equations formulated by Takahashi |S(J based on the string 
hypothesis. 

2. Non-linear Integral Equations (NLIE) found by Klu mpe r |pU and Destri-de Vega |}3| in the devel- 
opment of the quantum transfer matrix method |36|, |37fl, 

3. New integral equation derived by Takahashi 

The third one was recently discovered in an attempt to simplify the TBA equations pjj. It has also a close 
connection with the quantum transfer matrix |39] . The equation is explicitly given by 



u{x) = 2 + 
f = 



c \ x - y - 2i 

Tmu(O), 



exp 



2J/T 



(y + i) 2 + 1 



+ 



x - y + 2i 



exp 



2J/T 



1 dy 

it(y) 27ri' 



(3.3) 



where the contour C is a loop which counterclockwise encircles the origin. Numerically these three approaches 
provide perfectly the same data for free energy per site /. 
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Figure 2: P(n) at finite temperature 



4 Conclusion 

Let us mention again that the main result of this paper is the calculation of P(5) by means of the multi- 
integral representation. It is expressed only in terms of In 2 and Riemann zeta function with odd arguments 
with rational coefficients. This should be the general property of P(n). We expect that other correlation 
functions such as (SjSfe) will share this property. 

We have calculated numerically the value of P(n) and considered the asymptotic behavior. As was dis- 
cussed in ref. |nj, EFP P(n) shows a Gaussian decay at zero temperature, while it decays exponentially at 
finite temperature. 
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A Appendix 



Here we sho w the polynomials P , P x , P 2 which participate in the "canonical form" T5 given by the 
formula (2.5) 



( 5 ) _ 689 ,2,3,4 

"0 — -j^- A 2A 3 A4A 5 



(A.l) 
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(5) 796333 2 3299i a 13844 , 2 , a 165F7*. 3 . , 2 13463 2 , 3 ^ 

Pl = "15120" A4As _ ~6^ AlA4A5 + US - * 1 * 4 * 5 ~ "los"* 1 * 4 * 5 ~ ^512- A4A5+ 

1217491 , , 3 328189 i „ a 279917. 3 , 3 57619 , 2 . 3 543079 * . ,2,3 

^560~ AlA4A5 ~ ~M0~ AlA4As ~ ^45~ AlA4As + ^3T A4As ~ ^260~ AlA4A5 ~ 

9635, 2 , 2n 3 36199j . 3 . 2 . 3 811901 2 . 2 3 392107 2 3 125032 a 2 3 
— A!A 4 A 5 + ^^A 1 A 4 A 5 - 7560 A 3A 4 A 5 - ^ r A 1 A 3 A 4 A 5 + __A 1 A 3 A 4 A 5 + 

8.3, , 2 , 3 197459 4 39169 i . , 4 20873 2 4 66119* ,3, .4, 
7 A 1 A 3 A 4 A 5 - ~~^fgQ~ 4 5 - -T5lT XlX ^ ~ ^8T AlA4As + ^5^ AlA4A5+ 
66721i , 2 ,4 6301 , 2 4 5092H, 2 . 2 . 4 2045 3 a 4 1474 2 4 
"1512" 4 5 ~ ~90~ 1 4 5 + ~50T 1 4 5 ~ "8F AlA4A5 + — AsA 4 A 5- 

195492, , , 2 , 4 75367 ,2, ,2,4 , 21131* 3 ,2,4 , 10258 ,3,4 172369 i , 3 , 4 
— — — A1A3A4A5 ^-A!A3A 4 A 5 H — — A 1 A3M^5 H — 135~ A4 5 1512 - 

87209 2 3 4 , 21131* , 3 3 ,4 80497 i 3 4 12256 3 4 72071* a 3 4 

A!A 4 A 5 H Jgg - 1 4 5 504 3 4 5 27~ A l A 3A 4 A 5 H — AxA3A 4 A 5 + 

42262 ,3, , 3 , 4 161701 a sx4 ,29809* ,,2,3,4 ,29809 a a s 4 . . 

Ai A'AAaAk — A'MiAr. -h AiAyAiAr. -h Ai Aq A4. Ar. I /V.Z I 

189 5 1512 3 4 J 126 3 J 126 3 J v ; 

d(5) _ 41011 202861*. 3523 , , 6121 a 3433i, , 2 99 , 2 , 2 86029i, 3 

Fa " ~l^ + ^5^ A3 + ^ AlA3 + T2¥ A3 ~^^ AlA3 ~T AlA3 ~^890~ A3 ~ 

4813 , 3 130* 2 3 40,3,3 40799 2 , 596233 , , 5315 i . . 

AiAq -r Ai Ai -r Ai Aq — A5 — A3A5 ~\~ AlA3A5~r 

90 3 9 1 3 9 3 1890 5040 24 

582137i , 2, , 76147 , a 1640i , 2 , 2 , , 350 3 162292 , 3 2825 2 3 

— m , n A3A5 H r7-r-AiA 3 A5 — Aj^As H — A3A5 — — A1A3A5 — A1A3A5 + 

5040 120 3 9 36 4 

16602 ,3,3, 2993 ,2 , 1468283 2 , ,2 9213 a 35219 2 a 28679 i 2,2 
-^-A X A 3 A 5 - — A 5 + 504Q A 3 A 5 + -^-AiA 3 A 5 + ^-A 3 A 5 - _^AiA 3 A B - 

6715 2,2,2 151212 ,3,2 23695 3 a 3115* a 3 a 940 3 3 a , 47321* 3 

— ^— AiA 3 A 5 — A 3 A 5 — — AlA 3 A 5 H AxA 3 A 5 H — AiA 3 A 5 H ^260 

209881 , 3 9325 2 , 3 815092 , 2 , 3 15457 , a s 340* 2 . 2 . 3 4810 3 3 
T080 - " ^6~ AlA3A5 ~ ^04- A3A5 _ ^6~ AlA3A5 + "T" AlAaA5 " W A3A5+ 
1801* s 3 3055 3 3 3 3 3 3 6599 4 1151*. 4 53 4 1999 2 4 
— AiA 3 A 5 jg-A!A 3 A 5 + 160 2A!A3A 5 + ~^-X 5 ^j - 3 5 + 24~ AlA 3 A 5 _ "i~26 rA 3 A5_ 

599*, ,2,4 4225 ,2,2,4 113*. 3. 4 946 3 4 „.„., 2 , 3 ,4 , on ,3,3,4 



- A1A3A5 A? A 3 As — X3X5 g~ Al A 3 A s + 240 i\\ A3A5 + 80A? A^A 4 



(A.3) 
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